Sequential quantum-enhanced measurement with an atomic ensemble 
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We propose a quantum-enhanced iterative (with K steps) measurement scheme based on an 
ensemble of N two-level probes which asymptotically approaches the Heisenberg limit 5k oc 
jj-K/(^+i)^ i? the number of quantum resources. The protocol is inspired by Kitaev's phase esti- 
mation algorithm and involves only collective manipulation and measurement of the ensemble. The 
iterative procedure takes the shot-noise limited primary measurement with precision Si oc N'^'"^ to 
increasingly precise results 5k oc 7V~^'^. A straightforward implementation of the algorithm makes 
use of a two-component atomic cloud of Bosons in the precision measurement of a magnetic field. 

PACS numbers: 06.20.-f, 03.67.Ac, 03.65.Ta 



Estimating an unknown parameter of a quantum 
system usually involves a quantum probe prepared in a 
known state pq which, when brought into interaction with 
the system, evolves to a new state p^ under the action of 
a (so-called) quantum channeP Q^, po — )■ Q0(/5o) = P4>; 
by measuring a suitable observable of the probe in the 
state /5^, one can infer the value of 0. Due to the 
fundamental uncertainty of the quantum measurement, 
this information has a statistical character. In order to 
improve the knowledge on (p, one needs to repeat the 
measurement: measuring _R 3> 1 independent quantum 

probes pf^ -^ [Q0(/3o)]®'^ = /5^^ leads to n y/R in- 
crease in the measurement precision Sq — > So/yR — this 
is the standard quantum limit or shot-noise limit of mea- 
surement. Incidentally, nature provides us with a better, 
albeit ultimate^, Heisenberg limit 5q — > Sq/R, if one ex- 
ploits some of the quantum resources of the system under 
investigation. Exploiting such quantum effects enhanc- 
ing the measurement precision is the subject of quantum 
metrologjfflS. 

The quantum enhancement in the measurement preci- 
sion can be approached by using either parallel or se- 
quential strategies. In a parallel strategy, the origi- 
nal ensemble is divided into m sub-ensembles with n 
probes, R — nm, with each sub-ensemble prepared in a 
(maximally) entangled state [pq ] ™ -^ [Q<piPo )] *"; 
this results (ideally) in a ^/n enhancement of the pre- 
cision compared to the standard quantum limit, Sq — > 
So/^/mn = 6o\/rn/R, see Refs. El for the case of quan- 
tum interferometry and Refs. |6] for a general quantum 
channel. Alternatively, in a sequential strategy, instead 
of preparing entangled states, each of m separate probes 
is passed n times through the same quantum channel, 
pf™ — )■ [QS(Po)] , resulting in the same enhancement 
in precision, see Refs. W^. Roughly speaking, while parallel 
strategies make use of entanglement, sequential strategies 
exploit the coherent quantum dynamics as a resource in 
order to enhance measurement precision. While the se- 
quential strategies do not require creation of rather frag- 
ile entangled states, they do demand longer (coherence) 
times to allow for the completion of the measurement. 



Recently, quantum measurement protocols were im- 
plemented using two-component atomic ensembles with 
N particlegSffii'. In Ref. [SI interatomic interactions were 
used to create entangled, spin-squeezed stateJi^ of Bose- 
Einstein condensates (BECs) that allow one to go be- 
yond the standard quantum limit via a parallel strat- 
egy. Here, we suggest a sequential strategy without en- 
tanglement allowing to reach a given precision S within 
K - ln((5i/(5)/ln(\/iV) steps, where Si = Sq/VN is the 
uncertainty of the primary measurement and iV ^ 1 is 
the BEG atom number. This is achieved by subject- 
ing the BEG on each subsequent measurement step to a 
^ ViV-fold longer evolution under the unknown external 
field. Our approach resembles Kitaev's phase estimation 
algorithm^'^, however, it requires far less steps to com- 
plete the measurement due to our exploiting the large 
ensemble N of individual probes in the atom cloud. Our 
approach does not require any separate access to a par- 
ticular atom in the BEG and relies only on collective 
manipulations and measurements. 

Primary measurement. We wish to estimate the real 
angle (or phaseHl) g [— tt, tt] in an unknown unitary ro- 
tation f7z[0] = exp[— «(7z(/)/2] given an ensemble (with 
N ^ 1 probes) of spin-1/2 systems (qubits) and us- 
ing only collective unitary operations over the ensemble. 
We use Ramsey interferometry as the primary measure- 
ment and prepare all qubits in the ax — I polarized 
state (It) + ||))/V2. Applying Uz[(j)] to the ensemble 
results in the state [(e-"^/2|t) -h e+*'^/2|4,))/v^]®^. A 
Uy[~TT/2] pulse rotates the ensemble into the read-out 
state [cos (0/2) 1 1) + « sin((^/2) | i)]®^ and measuring the 
total polarization S^ = i^/N) J2i=i^^^ ^ o^^ arrives at 
one of the possible outcomes Sz — [N+ — N-]/N, with 
N+ e{0,...,N} and iV_ = N-N+ the number of qubits 
observed in the az — ±1 states. The probability to ob- 
serve a particular value Sz (we denote by X a realized 
value of the random variable X) is given by the Bernoulli 
distribution 
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conditioned on the unknown value of (p. Starting with 
an unbiased, homogeneous distribution P{4>) for the pa- 
rameter (j) (no a priori knowledge on (j)), the measure- 
ment of a particular value of Sz allows us to improve 
our statistical information on (p. Changing variables 
(j) ^f p = cos (f) and making use of Bayes' theorem, the a 
posteriori probability distribution P{p\Sz) of p S [~ljl] 
knowing the measured result Sz is enhanced by the factor 
V{Sz = Sz\(l))/'P{Sz), hence, after proper normalization, 



P{p\~S.) = 
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In the following, we consider the limit of large iV+ , N^ ^ 
1, where the distribution P{p\Sz) has a sharp peak near 
Sz] expanding the exponent in exp[A^+ ln(l-fp)-f A^_(l — 
p)] around the maximum, the distribution P{p\Sz) then 
can be replaced by the normal distribution 



P{p\~Sz)-- 



'2tt(j 



exp 



{p-Sz? 



2^2 



N 



(3) 



or in another notation, p ^ Af{Sz, cr^)- The overall result 
of the above ensemble measurement is summarized in the 
following statistical statement: Given the tolerance level 
/3 -C 1, the precision of the estimate p Ki Sz is given by 



PToh[\p~Sz\<gi(3)c7] =1-/3, 



(4) 



with g determined hy 1 — (3 — eii{g/\/2) and erf(a;) — 
(I/^/tt) /_ dte^* is the standard error function. Fi- 
nally, the distribution (|3| for p = cos (p provides us with 
the distribution function for the angle (j), 
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Complementary measurement. Since cos (j) is even in p, 
the estimate for (f> is distributed among two symmetric 
intervals 10 ± ^| < g/^fN with = | arccos(5z)|. Ex- 
panding ([5J) near (j) k, ztcj), the distribution of <j) is given 
by the sum of two normal distributions, 

^ ^ N{a4>,al), Ti = ^, (6) 
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describing two equiprobable alternatives a = ±1 for the 
angle (f) to be located near a4>. In order to distinguish be- 
tween the two alternatives of the primary measurement, 
we perform a second test by preparing the ensemble in the 

ay = +1 polarized state, ([|t) +«li)]/^/2) • Repeat- 
ing the Ramsey measurement, the second estimate for 
the parameter p' = sin should result in either -I- sin (j) or 
— sin 0, thus distinguishing the alternatives +(f> and —4> 
provided by the first measurement. Specifically, given 
the probability to observe a qubit in the cr^ = ±1 state 
conditioned on the result ap of the primary measurement 
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we find that the total polarization of the complementary 
ensemble is given by the sum of two normal distributions, 
S'^ ^ J2a-^(^'zaT'^'^)/'^' '^ith mean and variance 
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To construct an unbiased classification rule we define the 
regions E^ = {S-^l^i < S'z} and E+ = {S'z\S'z > S'z} 
with the boundary S'z set by the condition 



r{E^) 
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dS'zP{S'z)=V{E+) ^1/2. (9) 



In our symmetric situation, S'z = {S'z_ + S'z_^_)/2. Given 
a measured S'z we then assign the value a = I {a — 
— 1) whenever the event has been realized in E+ (E^). 
This assignment is prone to a misclassification error 
/3' = 'p(+\E_)ViE_) + P(-\E+)V{E+), where the con- 
ditional probability Vi+\E^) = V{E_\+)V{+)/V{E^) 
(and similar for V{—\E+)) follows from Bayes' theo- 
rem; for our unbiased classification rule, this reduces 
to V{a'\Ea) = V{Ea\a'). The conditional probabili- 
ties V{Ea\a') are easily obtained from the distributions 
r{S'z\a) = A/'(S';„, cr'2) and we find that 



/3' = [l-crf(|sin0|/%/2o-')]/2. 



(10) 



Away from the immediate vicinity of ~ (e.g., N+ > 5) 
and a typical number of probes N ~ 10'^, choosing 
N' ^ N results in a negligible probability f3' of misclas- 
sification. 

n-fold rotation. Next, we analyze the sequential (n > 
1-fold) application of the rotation Uz[<l>\ in the Ramsey 
measurement. Measuring the ensemble polarization Sz 
then provides an estimate for the parameter (/>„ = ncj). 
Given a measured result Szm the a posteriori distribution 
function for the parameter 0„ is given by Eq. ^ with 
Sz — ^ Szn and (p ^ (j)n, providing two Gaussian peaks 
at (/)„ = ±1 arccos(S'z„)| of width ex 1/^/N as a function 
of (/)„; the accompanying complementary measurement 
again selects one alternative a = ±1. 

A further estimate of the angle p is obtained by reading 
the distribution function Eq. ^ with cost/) — > cos(n0) 
as a function of 0. The periodicity of cos(n0) provides n 
different values of (j), all corresponding to the same value 
of cos(n0). The distribution function P„(0) for the angle 
(j) then has n peaks centered at 0„fc, 
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4>nk = —+2TT-, fc = 0,...,n-l, (12) 

n n 

each peak n times narrower than in the previous n — 
1 case. As a result, the n-fold measurement redis- 
tributes the original uncertainty 5i ~ cti among n dif- 
ferent equiprobable positions which we call the alterna- 
tives Ak. While this result does not give us more a 



posteriori information on the position of (j) than the 1- 
fold measurement (as confirmed by the Shannon entropy 
Hn = — / d(j)Pn(4>) lnP„((/)) coinciding for all n), the dif- 
ferent distribution of the probability allows us to gain in 
precision when combining the two measurements. 

Beyond the stiot noise limit. In order to take advantage 
of the n-fold measurement one has to identify the correct 
alternative among the n equiprobable distributions, see 
Eq. (11). This is done by combining the results of the 
1- and n-fold measurements. We define the interval /i — 
{4>\ \4> — 0i| < g{(3)ai} centered around the result (pi of 
the first measurement, see Fig. [T] with /3 the tolerance 
level of the first measurement, hence Prob[|(/> — (f>i\ < 
g{(3)ai] = 1 - /?, cf. Eq. Q. We call the alternative Ak 
compatible with the first measurement if the condition 

Pvoh[(j)eIi\(j)^M{4>nk,al)]>l~/3 (13) 

is satisfied. In order to satisfy the original confidence 
level in the second measurement, the maximum 0„fc be- 
longing to Ak must be located within the reduced interval 
/i< = {0||0-^i|<5(/3)(ai-a„)}. 



see 



Fig|T| 



By construction, the condition ( 13 1 is satisfied for at 
least one k, irrespective of the value of n. Choosing a 
small n, the gain in precision is small, hence we are in- 
terested in maximizing the value of n. On the other hand, 
for large n, the number of peaks compatible with Eq. ( 13 ) 
is larger than one and we cannot select the proper alter- 
native Ak. The optimal number riopt can be determined 
by considering the situation where the k-th peak is lo- 
cated at the left boundary of /j^, (pnk = (^i" 9{P){'^i~^ n) , 
while the next peak 4>nk+i is being pushed out from If 
across the right boundary of I^ with decreasing n, see 

Fig. [ij Obviously, when ^nk+i = 0i + gW){<^i - ct„) we 
still have two equally probable alternatives generating a 
large misclassification error /3 — 1/2. The task then is, to 
find the largest possible n compatible with a prescribed 
error /? ^ 1. 

After two measurements, a 1- and a n-fold, the a pos- 
teriori distribution function for the angle (j) is given by. 
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and Wk is the a posteriori proba- 



bility that the fc-th alternative has been realized, 

(^1 -<i>nkY' 
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For our arrangement (pnk = 4'i ~ giP)i'^i ~ '^n) and 
(l>nk+i = 01 - 9{P){cri - (Tn) + A„. Considering only 
these two peaks, the misclassification error of the k-ih 
alternative is given by /3 = Wk+i/{wk + Wk+i). Solv- 
ing for A„, we obtain the optimal number of rotations 
compatible with f3, n < riopt = l^iP, l3)/ai\ = n2 with 



;i + 21n[(l-/3)//3]/g2(/3)-i 

iy = TTg(p)- = — = . (16) 
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FIG. 1: Intervals Ji and If centered around 01. The peaks 
belonging to the alternatives Ak and Ak+i are shown for the 
extreme case where pnk coincides with the left boundary and 
pnk is on the verge of leaving If as n decreases, A,i — 2n/n. 



In the end, the precision estimate of two measurements 
with a 1-fold and a n2-fold rotation is given by 

Prob[|0-02|<5(/3)^2] = (l-/3)(l-/3), (17) 

where p2 is the measured and selected value of the pa- 
rameter p on the second step and the error bar is given 
by (72 = cTjia = o'i/n2. Thus the second measurement im- 
proves the precision by the large factor y/Nv{j3^j3) ^ 1, 
with i/(/3, /3) « 0.96 for /3 = /3 = 0.01 (for small /3, /3 < 1, 
vK (7r^2|ln/3|)/|ln;9|)[(l + ln/3/ln/3)i/2_i]. The in- 
fold iteration of this procedure will further improve the 
precision of the measurement. 

The above measurement protocol involves two sources 
of error, the estimation error /3 for the angle p to lie 
outside the interval li/) — 02| < 5(/3)f2, and the classifi- 
cation error /3 for an incorrect choice of the alternative 
Ak. While an estimation error can be ruled out on a 
subsequent step, this is not the case for a classification 
error. Indeed, for K = 2 and assuming that we have cor- 
rectly identified the peak p^k — p2i let us suppose that 
10 — p2\ > 5(/3)o'2, hence the true value p is outside the 
allowed range. Applying n^ — v j 02 rotations in the next 
step, none of the peaks pn^k will belong to the interval 
I2 defined through \pi — pi2\ ^ ff(/3)(<''2 ~ '^z) with prob- 
ability 1 — /3, signalling the error in one of the previous 
steps. The classification error cannot be caught in subse- 
quent steps: Consider two alternative angles p) localized 
either near a) the fc-th peak |0 — ^n^fel ^ 9{P)'^2 or b) the 
next peak \p - {(t>n2k + 27r/n2)| < ff(/3)o'2. Then apply- 
ing a ri3-fold rotation, the random parameter p^ = n^p 
will be localized within a) \p3 ~ n^pn2k\ < g{P)^j or b) 
103 — nzpn^k — 2'kv\/N\ < g{li)v. Since these intervals 
are largely overlapping, the hypotheses a) and b) cannot 
be distinguished. 

Iterating the process K times, our measurement pro- 
tocol satisfies the confidence criterion 



Prob[|0 -Pk\< g{P)<TK] = {\-P){i-P) 
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(18) 



with an — CTiici/j^]^^^ ^ N^^/"^. To reach a given 
precision 5 one then needs to perform Kg ^ 1 -|~ 
\n{aig / 5) / \n{u^/N) steps; at the same time, the over- 
all confidence level decreases exponentially oc exp(— /SX^) 



with the number Ks of steps. Moreover, an estimation 
error, i.e., none of the peaks belongs to the prescribed 
estimation interval, forces one to repeat the entire pro- 
cedure again. However, in practice VN ~ 30 — 40 and 
thus hardly more than a few steps are required to reach 
a good precision S. 

In counting the number of resources needed to reach 
the estimate (18), we choose as our basic unit the opera- 



tion Uz[(t)\ applied to a single qubit. The i-th step of the 
above procedure requires (for N qubits in the ensemble) 



R, = Nn, = N{v/a^y 



j\j-(i+i)/2 elementary 



operations. For v/ai S> 1, the entire K-step process 
then involves R = X]i;=i ^i ~ ^K operations. Express- 
ing the prescribed precision 6k {^ go'K = d'^iWi/'^Y'^^^) 
through the number R of resources, we arrive at the scal- 
ing 
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telling us that the Heisenberg limit is reached asymp- 
totically at large values of K. In an actual implemen- 
tation with ultracold atoms, the time needed to pre- 
pare the atomic ensemble is typically a few secondsSl, 
and increasing the number of rotations n does not sig- 
nificantly increase the overall duration of the measure- 
ment. The standard statistical measurement requires 
Kstd ~ {Sq/vNS)^ preparations to reach a precision S, 
which is exponentially larger than the Kg preparations 
required with the present protocol, Kstd ^ N^^~'^. 

Dephasing. In addition to the unitary rotation Uz [(j)] , 
the qubits may experience a stochastic field (p{t), e.g. due 
to uncontrolled interactions between qubits generating a 
different phase shift (p — >■ (j>+J dt ip{t) for each qubit. Av- 
eraging the singic-qubit density matrix over (p(t), the off 
diagonal amplitudes e exp{±i(f)) are reduced by the factor 
e — exp(— rTi/2) < 1, where we have assumed a Gaus- 
sian random field {(p{t)ip{t')) — rS{t — t') and ti is the 
exposure time of the primary measurement; the reduc- 
tion in these amplitudes after an n-fold rotation is given 
by e". The measurement of the ensemble polarization 
involves the parameter p — e" cos(n(/)) and the width in 
the a posteriori distribution function is cr„ — ?> ai/{ne^). 
The smallest attainable width aie ln(l/e) is reached after 
Uc = — 1/lne = Tc/ti steps, with Tc — 2/r the coherence 
time. 

Application. We analyze the use of our protocol to 
measure a constant magnetic field B with an atomic en- 
semble, considering a transition with differential mag- 
netic moment of order /is , the Bohr magneton. Assuming 
the prior knowledge that B < i3_|_, wc choose the interro- 



gation time of the first Ramsey sequence ti ^ 2nh/ ^^B^ 
such that the accumulated phase (j) = ix^Bti/H does not 
exceed 27r. This primary measurement results in a phase 
uncertainty of [5(l)\i = l/yiV, translating to a precision 
of [Ai3]i = hj \ji.-B,T\\[N in the field. In the following 
steps, the Ramsey time is increased as described above. 
The longest Ramsey sequence of duration t^ = ricTi pro- 
vides us with a precision of [(5(/)]niin = l/ricvN for the 
phase estimation, and thus a field precision of [ASJmin = 
h/fi^Tcy/N using K ^ 1 + ln(rc/Ti)/ ln(\/]V) steps. 

In a realistic situation the above procedure is feasible 
for small magnetic fields, since the duration ti of the first 
Ramsey sequence cannot be arbitrarily small — a typical 
Ti ^ 10~^ s corresponds to a field B+ ^ 1 G. With a typ- 
ical coherence time Tj, ~ 1 s and N ~ 1000 atoms one ar- 
rives after K = 5 steps at a precision AB ~ 3x 10~^ G. In 
order to measure higher fields one can exploit the phase 
periodicity and subtract an offset field. This requires 
prior knowledge that the field lies in an interval [B-,B+]. 
In this case, given a minimal time ti, we choose some 
field Bq S [i3_,i3+] that satisfies the matching relation 
lisBoTi/h = 2ttM with M the largest possible integer. 
The procedure described above is then used to measure 
the remaining small field b = B — Bq. 

Conclusion. It is interesting to compare our ensemble- 
based algorithm with Kitaev's original phase estima- 
tion algorithm involving individual qubits. In order to 
reach a prescribed precision 6, the latter necessitates 
K ^ \n{l/S)/ In 2 steps, a factor (In VN)/ In 2 larger than 
the ensemble-based protocol. For the Kitaev algorithm, 
the resources scale as i5 '^ InR/R (accounting for the 
fact that a final error probability f3 necessitates a smaller 
value P/K for the individual stepP^l) , which is better than 
our algebraic relation 6 ^ _R^^/^^+^^. Still, performing 
a few iterations on an ensemble of A^ ^ 10^ atoms with 
long coherence time appears as an attractive and practi- 
cal alternative to control isolated qubits over many steps. 

The sequential strategy discussed here is particularly 
useful in scanning probe measurements of a spatially 
varying field^. If the field distribution is a priori un- 
known, the duration of the initial Ramsey sequence has 
to be very short at each pixel of the image. Our scheme 
allows to optimally adapt the sequence at each pixel 
to reach fast reduction in the measurement uncertainty. 
This dramatically reduces the overall time to record a 
picture with a given precision. Finally, our ensemble- 
based sequential strategy could also be combined with a 
parallel strategy using squeezed BECs. 
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